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ABSTRACT 

In the era of synoptic surveys, the number of known gravitationally lensed 
quasars is set to increase by over an order of magnitude. These new discoveries 
will enable a move from single-quasar studies to investigations of statistical sam- 
ples, presenting new opportunities to test theoretical models for the structure of 
quasar accretion discs and broad emission line regions (BELRs). As one crucial 
step in preparing for this influx of new lensed systems, a large-scale exploration 
of microlensing convergence-shear parameter space is warranted, requiring the 
computation of O(IO^) high resolution magnification maps. Based on properties 
of known lensed quasars, and expectations from accretion disc/BELR modelling, 
we identify regions of convergence-shear parameter space, map sizes, smooth 
matter fractions, and pixel resolutions that should be covered. We describe how 
the computationally time-consuming task of producing ~ 290000 magnification 
maps with sufficient resolution (10000^ pixel/map) to probe scales from the in- 
ner edge of the accretion disc to the BELR can be achieved in ~ 400 days on 
a 100 teraflop/s high performance computing facility, where the processing per- 
formance is achieved with graphics processing units. We illustrate a use-case for 
the parameter survey by investigating the effects of varying the lens macro-model 
on accretion disc constraints in the lensed quasar Q2237+0305. We find that al- 
though all constraints are consistent within their current error bars, models with 
more densely packed microlenses tend to predict shallower accretion disc radial 
temperature profiles. With a large parameter survey such as the one described 
here, such systematics on microlensing measurements could be fully explored. 

*Research undertaken as part of the Commonwealth Cosmology Initiative (CCI: www.thecci.org), an 
international collaboration supported by the Australian Research Council. 
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Subject headings: gravitational lensing: micro - quasars: individual (Q2237+0305) 



Introduction 



Cosmological gravitational microlensing refers to the regime of gravitational lensing 
where individual stars within a lensing galaxy act to magnify a background source, usually 
a quasar. The most obvious effect of microlensing is to produce an uncorrelated bright- 
ness change with time in a single image of a multiply-imaged source; an intrinsic vari- 
ation in source flux would appear in all images separated by the individual image time 
delays. Cosmolog ical microlensing was first observed in the lensed quasar Q2237+0305 by 



Irwin et al.l ( 119891 ). and it is likely to be taking place on some l evel in all qu a sars m ulti 



ply imaged by a foreground g alaxy. 
Schmidt fc Wambsgansj ( 120101 ). 



For recent summaries, see IWambsgansd (120061 ) and 



While gravitational microlensing is achromatic, the magnitude of microlensing-induced 
fluctuations depends strongly on the size of the emission region. Small sources (relative 
to the Einstein radius of the microlenses) are affected more strongly than larger sources. 
Detections of microlensing signals in a single waveband can therefore be used to constrain 
the size of the quasar emission region. For example, optical obse rvations of Q2237-I-0305 
sugges t a V-band accretion disc radius of Ry = G-'^-l'^ x lO^^cm (jPoindexter fc Kochanek 

mm- 



Multi-wavelength observations present an opportunity to constrain the wavelength- 
dependent structure of quasar emission regions. First constraints on the temperature profile 



of quasar accretion dis c s have already been pub l ished (lAnguita et a. 



Eieenbrod et al. 


2008 


Mediavilla et al. 


2011 



Poindexter et al. 


20C 


8 


Mosquera et al. 


2011a 


)• 



Floyd et al 



120091: 



120081: 



Bate et al. 



Blackburne et al 



2008 



2011 



Interpretation of broad emission l ine region (BELR] sign als are more comple x, however 
some early attempts have been made. iLewis et al.l ( 119981 ) and IWayth et al.l ( 12005! ) measured 
the relative sizes of the CIII] / Mgll BELR and cont inuum emission region in Q2237+0305 us- 
ing microlensing observations. I Aba j as et al.l (120071 ) proposed a biconical BELR as a possible 
explanation for a recurrent blue- wing enhancemen t in the emission l ines o f SDSS J1004-I-4112 
( iRichards et al.l 12004 : iGomez- Alvarez et al.l |2006| ) . lO'Dowd et al.l ( 12010[ ) compared a differ- 
ential microlensing sign al observed across the CIII] emission line in Q22374-0305 with various 
BELR geometries, and ISluse et al.l ( 120111 ) studied the broad emission lines in 39 epochs of 
Q2237+0305 data obtained with the VLT. 
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While ~ 100 multiply imaged quasars have now been discov eredjll only ~ 20 of these 
have measured light travel time delays between lensed images - see lOguril (120071 ) for a recent 
compilatio n. New, high-cade nce, all- sky surveys such as those planned for the SkyMapper 



Teles cope (IKeller et al.l 120071 ) and the Large Synoptic Survey Telescope (LSST; llvezic et al. 



20081) are likely to disco ver a large number of gravitationally lensed quasars. For example, 
Oguri fc Marshall! (|2010[ ) predict that LSST will find approximately 8000 new lensed quasars 
across its 10- year 20000-deg^ survey, 3000 of which will have good measurements of the time 
delays between images. 

The anticipated increase in the number of known lensed quasars will take measurements 
of quasar structure beyond the realm of single-object studies and into that of statistically 
meaningful samples. For example, we can test whether a single accretion mechanism is 
consistent with all obs ervations, exa r aine s caling relations between black hole mass and 
accretion disc radii [see [Morgan et al.l (120101 ) for such an analysis with a sar nple of 11 lensed 



quasar s], and investigate secondary parameters such as Eddington ratio (e.g lMosquera et al. 
201 Ibh . 



Analysis and interpretation of this wealth of new information regarding quasar emission 
regions will require extensive simulations. The basic tool for microlensing is the magnifica- 
tion map: a pixellated approximation to a region of the (background) source plane, with a 
magnification value defined at each physical source location. By convolving magnification 
maps with wavelength-dependent geometric model profiles, both simulated lightcurves and 
statistical magnification distributions can be extracted, ready for comparison with observa- 
tional data. 

The characteristic caustic structure of a magnification map depends on the convergence 
and shear (described in Section [2]). These parameters are local properties of the global macro- 
lens model at the position of each lensed image. Since each lensing environment is unique, 
otentially ^ 6800 new combinations of convergence and shear will need to be simulated 
Oguri fc Marshall! ( 120101 ) predict ~ 14 per cent of ~ 3000 LSST-discovered systems will be 
quadruply imaged, with the remaining systems dominated by double images, resulting in a 
mean image number of 2 .28]. Moreover, macro-lens models are not uniquely determined. 



For example. IWitt et al. ( 119951 ) modeled lensing systems with three different lens potentials 
and Wyithe et al.l (2002) contains a compilation of lens models for Q2237+0305 by eight 
different authors: each of which predicts different convergence and shear values for the 
lensed images. The potential need for multiple models, and hence mangification maps, for 
each newly discovered system further increases the expected coverage of convergence-shear 



^See the CASTLES Survey webpage: | http : / / www . cf a . harvard . edu/glensdata/| 
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parameter space. 



Although th e computational time needed fo r analysis of magnification maps can be 



considerable (e.g. iPoindexter fc Kochanekll2010al ). generating sufficient high-resolution, sta- 



tistically useful magnification maps is itself a computationally-demanding task. Over the last 
two decades, the most widely-used approaches to generating magnification maps are based 
on backwards ray-shooting (described in more detail in Section 2). Here, light rays are prop- 
agated from the observer, through the lens plane, requiring the calculation of defiections by 
individual lenses, and mapped onto the pixellated source plane. This computation is 
repeated for a sufficiently large number of light rays, typically O(IO^) or higher, in order to 
achieve statistically significant coverage per pixel in the source plane. 



In the hierarchical tree method (jWambsgansslll990l : IWambsgansslll999l ) . individual lenses 
are approximated by higher-mass pseudo-lenses, depending on their distance fro m a given 
l ight r ay, thus avoiding the requirement for iV* defiection calculations per light ray. iKochanek 
(I2OO4J ) describes a computationally-efficient approach to reducing the effect, where the 
lens and source planes are treated as spatially periodic, allowing the use of Fourier meth- 
ods (similar to the P^M me thods used for niany-p article gravitational force computations). 
Mediavilla et al.l (120061 ) and lMediavilla et al.l (l201l[ ) adaptively tessellate the lens plane, with 
the lens equation used to map the lattice of polygons from the image to source plane, where 
the polygons are efficiently mapped onto source pixels. This approach reduces the total 
number of light rays required to reach a given statistical accuracy by over three orders of 
magnitude. 

The emergence of the graphics processing unit (GPU) as a means of accessing low 
cost, high performance, massively parallel computing, presents an alternative opportunity 
for speeding-up defiection angle calculations with an algorithmically-simple solution. As- 
tronomers have been amongst the early adopters of GPUs for scientific computation, report- 
ing typical speed-ups between x30 to xlOO over what they can achieve on single (or few) 
core CPUs. While further wo rk is required to fully unde rstand which types of computations 
are best suited for GPUs [see iBarsdell et al.l (I2OIOI ) and iFluke et al.l (1201 ll ) for a discussion 
of these issues specifically aimed at astronomers], algorithms that exhibit the properties of 
high arithmetic intensity and a high-level of data parallelism are the most promising targets. 



Thompson et al.l (120101 ) demonstrated that by taking advantage of the highly parallel 



GPU- architecture, a "brute force" implementation of ray-shooting is now a practical alterna- 
tive. Moreover, Bate et al. (2010) have shown that on current generation multi-GPU systems 
(e.g. the 4-GPU NVIDIA Tesla S1070 unit with a processing peak of ~ 2.5 terafiop/s for 
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single precision computatioiiE), the direct code on GPU executes in runtimes comparable 
too, and often faster than, the single-core tree code. 

Our target system is the Australian GPU Supercomputer for Theoretical Astrophyiscs 
Research (gSTAR§), a hybrid CPU+GPU cluster containing 102 x NVIDIA C2070 GPUs 
and 21 x M2090 GPUs, connected via QDR Infiniband, with a (theoretical, single-precision) 
peak of ~ 130 teraflop/s. Throughout this work, we estimate run-times for two different 
GPU cluster configurations: we assume a slightly more conservative 100 teraflop/s mode of 
operation for the phase one gSTAR system, gSTARioo, and a somewhat optimistic future 
upgrade with a peak of 600 terafiop/s, gSTARgoo, based on the original system design goals. 
In both cases, we expect to achieve at lea st 50 per cent of the theoretical peak, based on the 



benchmarking in [Thompson et al.l (120101 ) . 



Our goal in this paper is to demonstrate that a high resolution, large-scale, microlensing 
parameter survey is feasible on a GPU-enabled computing cluster. This requires generating 
statistically useful samples of high-resolution magnification maps over the entire range of 
convergence-shear parameter space - a theoretical equivalent of an "all-sky" survey. By 
providing ~ 290000 pre-computed magnification maps as a Web-accessible resource, rapid 
first look computations would be possible as new microlensed quasars are discovered. High- 
resolution coverage of parameter space enables investigation of degeneracies between macro- 
model parameters and inferred plausible models for quasar emission regions, and presents 
new opportunities to test (e.g. via mock observations) theoretical models of the accretion 
disc and BELR well before the new era of microlensed quasar discoveries begins. A parameter 
survey does not remove the need for detailed modelling of individual systems: for example, 
simultaneous multi-image fits to light curves including dynamical motions of the stellar lens 
populations, requiring g e neratio n of magnification maps for each observation epoch, as per 



Poindexter fc Kochaneki (l2010al l 



This paper is organised as follows. In Section [2l we discuss the theoretical consider- 
ations for a survey of the entirety of convergence-shear parameter space on a GPU-based 
computing cluster. We discuss magnification ma p size and resolut i on con siderations, and use 
simulation time estimates based on the results of [Thompson et al.l (120101 ) to inform a realistic 
simulation strategy. In Section E] we illustrate one possible use for our survey by examining 
the systematic effects on accretion disc constraints in the lensed quasar Q2237+0305 when 
the lens macro-model is varied. We present our conclusions in Section HI 



^1 teraflop/s = 10^^ floating point operations per second. 

■^The gSTAR facility was part-funded by Astronomy Australia Limited through the Australian Federal 
Government's Education Investment Fund. 
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Throughout this paper, a cosmology with Hq = 70kms ^Mpc ^, flm = 0.3 and Q\ = 0.7 
is assumed. 



2. Survey design 



In this section, we discuss the theoretical considerations that inform the choices of 
parameters for a high-resolution, cosmological microlensing parameter space survey. We 
discuss the ways in which observational constraints shape the magnification map sizes and 
resolutions that are required to produce a useful resource for the microlensing community. 

There are two techniques for constraining qua sar structure from microlensing obse rva- 
tions: analysis of time-dep endent light curves fe.g. IWyithe et al.ll2000t iKochanekl l2004j ) . or 
single epoch observations ( IWitt et al.lll995l : iBate et al.ll2008l ). Both rely upon a lens model, 
which is used to generate magnification maps that describe the fluctuation of magnification 
induced by the microlenses, projected onto the source plane. The lens model is defined in 
terms of the (external) shear, 7, and the convergence, k. Shear describes the distortion 
applied to images and is assumed to be dominated by matter in the host galaxy, however 
any matter along the line of sight to the source can contribute. Convergence describes the 
focussing power of the (macro-) lens, and includes contributions from both smooth matter. 
Kg, and compact objects, according to k = Ks + k^. We define the smooth matter fraction 
to be 

S = Ks/k. (1) 



Once a set of model parameters have been defined, a magnification map can be generated 
using the gravitational lens equation in the following form: 



1-7 

1 + 7 



K^X - ^ 



i=l 



(X - Xj) 



X — X,- 



(2) 



which relates the two-dimensional location of a light ray, x, with a source position, y, for 



7l{M) 



(3) 



compact lenses distributed in an ellipse of area A, at lens plane coordinates, x,, and with 
mean lens mass, (M). 

Using Equation ([2]) to shoot light rays at known positions backwards through the lens 
plane to the source plane, we can count the number of rays reaching each pixel, Nij. Com- 
paring these values to the average number of light rays per pixel if there were no lensing. 
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Nrays, we obtain estimates for the per pixel magnification: 



(4) 



For large A^* and large total number of light rays, the direct computat ion of Equa- 



tion (1 


2]) on a sing] 


e-core 


1990 




Wambseanss 


1999) 



The hierachical tree-code (jWambsganss 



Wambsgansd Il999l ) overcomes this limitation by replacing the A^* lenses with a re- 
duced n umber of higher -mass pseudo-lenses. A similar approach is used in the Fourier 
method (jKochanekl |2004J ) . where the lens and source planes are treated in a periodic fash- 
ion, which helps reduce the total number of light rays required, and also removes the edge 
effects, making raore of t he magnification map u sable for analysis. The tessellation approach 
(iMediavilla et al.l l2006l : iMediavilla et al.l 1201 ll ) reduces the total number of light rays by 
using a sophisticated mapping of polygonal regions from the lens plane to the source plane. 

While all of these algorithms lead to improved computational efficiency, they do this at 
the expense of algorithmic (and hence implementation) simplicity. The high level of paral- 
lelism in the direct ray-shooting algorithm, coupled with the arithmetic intensity inherent 
in Equation ([2]), are ideal attributes for implementation on a GPU. For details of the GPU- 
D code, including optimisation issues and timing tests, see [Thompson et al.l (120101 ). For 
timing and acc uracy comp a risons between GPU, hierarchi cal tree codes and parallel large 
data codes, see iBate et al.l (120101 ). IMediavilla et al.l (120111 ) report on a timing comparison 
betw een an imp r oved implementation of the tessellation method and the GPU-D results 
from lBate et al.l (120101 ): they find a two order of magnitude improvement in computational 
efficiency, but the total run-times of the two approaches were comparable. 

As the three (currently) CPU-only approaches for generating magnification maps share 
a high-degree of parallelism, particularly through the need to calculate defiections of many 
light rays, opportunities exist for their implementation on massively parallel architectures. 
Such work is beyond the scope of this present paper. We adopt the viewpoint that we have 
a working GPU ray-shooting code, and have access to a ~ 100 terafiop/s supercomputing 
facility, gSTAR, and hence present what could be achieved now. Any progress in reducing 
the total computation time should be traded off with t he additional amo unt of time required 
to implement an alternative, more complex, solution (jpiuke et al. 2011 ). 



2.1. Map size and resolution 

In microlensing simulations, we run into the usual conflict between map size (measured 
in units of the Einstein radius) and map resolution. We need maps that have sufficient 
resolution to resolve the very inner regions of quasar accretion discs (~ lO^^cm; see for 
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example the analysis of RXJ1131-1231 in iDai et al.l |2010[ ). and sufficient size to map out 
quasar broad emissio n lines (~ lO^^cm fo r the CIII]/CIV/MgII lines; see for exa naple the mi- 



crolensing analysis of lWayth et al.ll2005l . or the reverberation mapping analysis of lKaspi et al. 



20071 ). One option is to provide two sets of maps with sizes and resolutions tuned to the 
extreme physical scales for the accretion disk a nd BELR. In light of recent multi- wavelength 



studies of accretion di s k sizes and structures (jPoindexter et al.l l2008l : [Morgan et al.l 12010 



Blackburne et al.l 1201 ll : iMosquera et al.ll2011al ). the approach we adopt is a single map that 
spans both physical scales. This permits the use of self-consistent quasar models for studying 
future, simultaneous multi-wavelength observations extending from accretion disk to BELR 
scales. 

To set the pixel resolution required, we take the last stable orbit of a Schwarzschild 
black hole of mass, Mf,h: 



R 



6GM, 



bh 



ISCO 



8.86 X 10^^ 



M, 



bh 



V lO^M, 







cm, 



(5) 



where G is the gravitational constant and c is the speed of light. We define this value for a 
IO^Mq black hole as the fiducial radius, Rj^sco- ^ fixed map size, the pixels required to 
reach exactly this resolution will vary from system to system, as the Einstein radius depends 
on the angular diameter distances, -Djj, between the observer, o, lens, and source, s. 
Projected onto the source plane, the Einstein radius is: 



^0 



od 



(6) 



for a microlens with mean mass, (M). For a background quasar at redshift z = 2.0, a lens 
at z = 0.5, and a mean microlens mass (M) = O.SMq, the Einstein Radius is 2.87 x lO^^cm. 

How, then, do we determine the required number of pixels? We start by considering 59 
of the currently known lensing systems with redshifts available for both sour ce and lens using 
data f rom the CASTLES Survey (this is a more conservative choice than iMosquera et al. 
2011bl . who used estimates of lens redshifts to build up a sample consisting of 87 systems). 
Choosing (M) = O.SM©, we calculate the mean Einstein radius (r^o) and its standard devia- 
tion to be: 



{r]o) = (2.93 ± 1.20) x W^h 



-1/2 



70 



(M) 



1/2 



cm. 



(7) 



.O.SMq^ 

The Einstein radius of each system, along with the mean and standard deviation, are plotted 
in Figure dJ To achieve our target pixel resolution, R{^sco^ maps would require an average 
of ~ 330 pixels per rjQ. 
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The two systems with unusually large Einstein radii are Q2237+0305 and MG1549+3047. 
In both cases, the large Einstein radii are a result of lensing galaxies at unusually low red- 
shifts: Zd = 0.04 for Q2237+0305, and z h. = 0.11 for MG 1549+3047. The latter is also 
an extended source seen as a radio ring (ILehar et al.lll993l ). and is not known to exhibit 
micro lensing, so it is an atypical object. Removing these two systems essentially does not 
change the measured mean Einstein radius, but reduces the standard deviation by a factor 
of 1.79. We make no comment here on selection effects that might bias this measurement; we 
simply note that the known sources are a reasonable place to start. Existi ng gravitationall y 
lensed quasars have been discovered using a wide variety of techniques (see lKochanekil2006l ). 
which are likely to be developed and refined in the era of large synoptic surveys. 

Next, we determine the preferred physical size for the magnification maps. There are 
two main considerations here: that they be large enough to conduct simulations of quasar 
BELRs, and that they are wide enough to allow long light curves to be simulated. 



Assuming the quasar H/3 BELR radius-luminosity relati on (iBentz et al.ll2009l) . Rbelr oc 
2^0.519^ holds for high- luminosity, hi gh-redshift quas a rs (see iKaspi et al.l 120071 ) . we can esti- 
mate the expected BELR radius. iMosquera et al.l ( l2011bl ) performed this calculation for 
their sample of 87 lensed quasars; for the 44 that overlap with our sample, the estimated 
mean H/3 BELR and its standard deviation are (1.71 ± 1.47) x 10^'' cm, or (5.84 ± 5.01) (r^o). 

Under the usual criterion that sources with sizes smaller than or approximately equal 
to the Einstein radius can be significantly microlensed, this suggests that H/3 BELRs are 
generally too large to undergo significant microlensing. However, the BELR is known to be 
stratified, with higher ionisation lines emitted closer to the central source. In NGC 5548, for 
example, th e CIV line is ernitted from a region a fact or of two smaller than the H/3 line (see 
for example iPeterson et al.lll99ll : iDietrich et al.lll993l ). Microlensing estimates of the size of 
the GUI] emission region in the lensed quasar Q2237-I-0305 suggest a radius similar to th e 



Einstein radius for that source (IWavth et al. 



2005 



O'Dowd et al. 



2010 



Sluse et al.ll201lh . 



In line with the microlensing analyses of lAbajas et al.l (120021 ) and iLewis &: Ibatal (|2004[ ) , we 
therefore assume an approximate outer radius for the high ionisation BELR of Ir/o- 

A map with side length bOrjQ, convolved with a Ir/o radius BELR, would provide ~ 500 
statistically independent data points [taking int o account edge e ffects due to convolution - 
these can be avoided using the Fourier method ( iKochanekl |2004| )] . This drops dramatically 
if the BELR is much larger than Irjo, however repeat generation of statistically independent 
maps allows us to compensate for this drop off. A 50^7^0 map would require ~ 16500^ pixels 
in order to resolve Rj^sco- 



Fortunately, it appears that the effective velocities of lensed quasars are likely to be 
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small enough that a map large enough to enable simulations of qu asar BELRs will also be 
large enough to allow adequate sampling of microlensing lightcurves. iMosquera et al.l (l2011bf ) 
estimated that the average Einstein radius crossing time (the time taken for a source to travel 
a single Einstein radius) for a sample of 90 gravitationally lensed quasars peaks at ~ 23 years, 
with a range of 8 to 44 years, assuming a mean microlens mass of (M) = O.SMq. We can 
therefore simulate observations of duration ~ 10 years comfortably with a magnification map 
large enough for BELR simulations. 



2.2. Timing 

To further refine our choice of p arameters, we inv e stigat e anticipated run times. Using 



the empirical relationship, Tqpu, in [Thompson et al.l (120101 ) for a 4-GPU NVIDIA S1070 
Tesla unit: 



Tgpu = 1.4 X loM ) ( ^ ) ( ^ I sec 



^40962 J \10y V 100 

we can estimate the likely run time for any given magnification map. We have not estimated 
times for t he hierarchi c al cod e, but note that we expect them to be longer overall based on 



Figure 3 of Bate et al.l fl2010h . 



Thompson et al.l (l2010[ l obtained a sustained processing performance of 1.28 terafiop/s 
on the S1070, which is approximately half the nominal peak performance (2.488 terafiop/s). 
We note that quoted peaks are rarely achieved in scientific computations, as they typically 
require more than a dual-issued multiply and addition per GPU clock cycle. Assuming 
this half-peak performance can also be obtained with gSTARioo (gSTARgoo), we can expect 
approximately 50 (300) terafiop/s performance utilising the entire facility - an optimistic 
viewpoint that ignores real-world scheduling issues and the needs of other gSTAR users. 

Following Bate et al. (2010), we set Arrays = 1000, as this was found to be sufficient for 
the parameters explored in that paper, and A'pix is allowed to vary in line with the discussion 
in Section [TTl The remaining unknown in Equation (|H]) is the number of microlenses, A^,,, 
for each (k, 7) pair, which we calculate from Equation 

We also need to determine the total number of (k, 7) combinations to cover parameter 
space. To avoid making any assumptions about the particular lens system each parameter 
combination represents, we suggest that unifom tiling of k-j parameter space is the most 
conceptually appealing approach. It is likely that a single {k, 7) combination will be relevant 
for more than one lensing system due to degeneracies in rjo. By examining published macro- 
model parameters for known microlensed quasars (see Table [1] and red symbols in Figure [2]) 
we see that choosing a maximum k = 7 = 1.7 gives us good coverage. We choose uniform 
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increments Ak = A7 = 0.01, consistent with the accuracy of most of the quoted macrolens 
models. This results in a 170 x 171 sampling of K-7 parameter space {k = 0.0 is excluded, 
as in this case there are no microlenses) . 

We obtain timing results for 10 values of the smooth matter fraction s = Kg/ 1^ per (k, 7) 
combination. The smooth matter fraction ranges from < s < 1 in increments of As = 0.1. 
For a given (k, 7) combination, increasing the smooth matter fraction decreases the number 
of stars A^^,. This is accounted for in our timing results. 

We note that although we have attempted to collate all the (k, 7) combinations in the 
literature for Table [H the list may not be complete. In addition, many authors publish 
enough information on lens macromodels to extract c onvergences and shears , but do not 



quote the values explicitly (see for example Table 6 in iBlackburne et al.ll201l[ ). We do not 
include such models in our analysis. 

Timing contour plots are provided in Figures [2] and |3l The grayscale in both figures 
ranges from hours (white) to greater than or equal to 24 hours (black) . Note that these are 
actual hours, not device hoursjj Overplotted in both figures (red stars) are the models from 
the literature (see Tabled]). Figure [2] shows the complete timing; the white line indicating 
points where (1 — k — 7) or (1 — k + 7) equal zero corresponds to magnification maps that 
are formally infinite in size, and hence are excluded from our computation. Figure |3] shows 
timing coverage with regions requiring > 1 gSTARioo day to process masked out. 

In Table [21 we present six possible strategies for a GPU-enabled microlensing parameter 
survey. The first, labelled 'Complete', assumes that there are no limitations on the time we 
can run the survey. Parameter space is completely covered, at the map size and resolution 
that best suits the large microlensing simulations we would ideally like to run (see Sec- 
tion Obviously, this parameter set is impractical; it would require running gSTARioo 
(gSTARgoo) non-stop for ~ 115 (~ 19) years! 

In the 'Realistic' strategy, we sacrifice both map size and pixel resolution in order to 
speed up the time to complete the simulations. In addition, we exclude from this strategy 
any magnification map that requires more than 1 gSTARioo (gSTARgoo) day to simulate. 
This amounts to 739 (215) excluded combinations of k and 7, or ~ 4 (~ 1) percent of the 
total number of simulations. The simulation time is a more reasonable ~ 414 gSTARioo days 
(~ 167gSTAR6oo days). 

Simulation strategy 'Stage 1' gives an indication of the time required to produce one 



*These are based on our estimate of processing on the entire cluster at 50 (300) teraflop/s, so that 1 
gSTARioo (gSTARgoo) hour corresponds to ~ 10^^ 10^^) computations in one real hour. 
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magnification map for eacli of tlie model parameters in Table [H All of the models in that 
table, representing most (if not all) microlensing models in the literature, would be simulated 
in ~ 122 gSTARioo days (~ 40 gSTARgoo days) using GPU-D. Starting with this strategy, 
successive processing stages can then fill in the rest of parameter space, moving outwards in 
regions centred around the current known models (Stages 2-4). 

We note that a significant amount of simulation time is spent in maps where (1 — k — 7) 
or (1 — K + 7) are close to zero. We have chosen not to include maps with simulation times 
greater than 1 gSTAR day in our timing. If this cutoff were reduced to 1 gSTAR hour, the 
'Realistic' strategy would take ~ 80 gSTARioo days (~ 34 gSTARgoo days), with 1617 (657) 
fewer {k, 7) combinations. 



2.3. Data products 

Considering the total (unmasked) K-7 parameter space sampled 170 x 171 = 29070 
times for a single smooth matter value, and with no repeat maps, the total uncompressed 
data volume is ~ 53 TB. With consideration of appropriate data types, file formats, the 
use of compression, and an appropriate database infrastructure to handle queries and map 
retrieval, we assert that data storage and service is reasonable. 

Magnification maps do not have negative values, so we can use an unsigned data type. 
For the parameters we propose, an unsigned short integer (16-bit) should be sufficient in most 
cases (for typical A^'rays ~ 1000), as this gives us a range from 0. . . S'max = 65536 discrete 
values. If there are cases where individual pixel values are Nij > S'max, it is unlikely that there 
are S'max distinct magnification values. In such a case, pixel values can be remapped into the 
range . . . S'max, and an index array can be built and stored as metadata. In uncompressed 
form, and writing maps in a binary file format, a single A^pix = 10000^ map requires 191 MB 
of storage]^ 

To examine the likely effects of compression, we generate ten maps with randomly as- 
signed values (up to S'max) for each pixel. We then trial two standard Unix compression tools, 
gzip and bzip2, and record the compression ratio and the times to compress and decompress. 
We find that gzip compresses the random maps consistently to 171 MB with compression 
and decompression times below 16 sec and 5 sec respectivelyj^ Compare this with 157 MB 



^We use computing notation where 1 MB = 1 megabyte = 1024 kilobytes, and similarly for gigabyte 
(GB), terabyte (TB) and petabyte (PB). 

^Measured on a 2.8 GHz, Intel Core 2 Duo, MacBook Pro, so these are upper limits. 
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file sizes for bzip2, and < 38 sec for compression and < 17 sec for decompression. 

We expect compression of real maps to be somewhat better than random pixels, as 
the range of per pixel values should be lower than a uniform distribution over C.S'max- The 
presence of caustic structures means that there are sequences of repeated values in the maps, 
which tends to aid with the run- length encoding of gzip, and the block-based compression 
of bzip2. We tested compression of a magnification map with k = 0.416, 7 = 0.471, 



corresponding to the /3 = 1.0 model for image A in Q2237+0305 in lWitt et al.l (119951 ). For 
a 10000^-pixel map, and Arrays = 1000, the count range was 15 < Nij < 7060. In this case, 
gzip and bzip2 reduced the file size to 113 MB and 80 MB respectively. Times to compress 
and uncompress the map were 23 sec and 5 sec for gzip and 26 sec and 13 sec for bzip2. 
Despite the slightly longer decompression times, which is actually the more critical factor 
for accessing data stored on disk in compressed form, our testing suggests that bzip2 is the 
preferred alternative, with likely compression ratios between 40 per cent and 78 per cent. 
These estimates are consistent wi th the compression ratios of 2.5 (40 per cent) to 3 (33 per 



cent) for OGLE data reported by iPoindexter fc Kochanekl (j2010al ). 



The overhead in decompressing any given map for analysis is small compared to the 
time taken to generate the map or to transfer large numbers of maps over Ethernet. As the 
cost per TB of data storage continues to decrease at a rapid rate, we conclude that data 
storage and access, even with repeat maps and sampling of smooth matter values, is not a 
significant concern. An efficient. Web-based data delivery mechanism would be required to 
serve sets of maps on the basis of user requests, and we are currently considering options. 



3. First Application 



We now present an example of one possible use for a microlensing parameter survey. 
In this section we analyse the effect of varying lensing parameters (specifically k and 7) on 
accretion disc constraints. 

As discussed in Section 1, lens macro- models are not uniquely determined. In systems 
where a wealth of observational data is available, complex models of the mass dis tribution in 



the le nsing galaxy can be constructed (see for example the Q2237+0305 model of iTrott et al. 



20101). More often, relatively simple parametric models are used to describe the lens mass 
distribution. These models are typically found to fit the lensed image positions adequately 
within observational errors, although a dditional mass components (perhaps due to other 
nearby galaxies, e.g. iPooley et al.l 120061 ) are sometimes required. 



The main uncertainty in these models is the percentage of the surface density in smoothly 
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distributed matter. This uncertainty has been included in many microlensing analyses, 
either as a free parameter in accretion disc measure ments, or as the paranaeter of interest in 



analyses of lens da rk matter fractions (see especially iMediavilla et al.ll2009l : iBate et al.ll2011 
Poolev et aliboilh . 



This uncertainty is typically confronted in one of two different ways. In the first tech- 



nique, the galaxy is mode 
White (NFW) halo (e.g. 


led with a de Vaucoulers profile embedded 


in a Navarro-! 


^'renk- 


Morgan et al. 


2006 




Poindexter et al. 


2008 




Chartas et al. 


2009; 


Morgan et al. 


2010; 


Dai et al. 


2010 


). The mass 


to light ratio is varied, typically from 1 (no 



NFW halo) to 0.1 (NFW halo dominates the mass distribution). This results in a sequence 
of lens macro-models, from which microlensing convergences and shears can be extracted. 

In the second technique, a single macro-model is produce d (usually a singular isothermal 
ellipse with external shear) to describe the lensing galaxy (e.g. Bate et aPboOS ; Mediavilla et al 
20091 : iBlackburne et al.l 1201 ll ) . A single set of microlensing convergences and shears are ob- 
tained from this model. The unknown smooth matter fraction is then accounted for by 
allowing the smooth and stellar components of the convergence to vary, subject to the 
constraint k = + n^. We have assumed this second technique in our discussion on survey 
design, however we note that the magnification maps generated by such a survey are equally 
useful regardless of lens modelling technique. 

Since a variety of lens modelling techniques are used, each providing acceptable fits 
to lensed image positions, it is reasonable to ask whether source quasar parameters can 
be uniquely determined from microlensing observations. We present a first attempt at a 
comparison between accretion disc constraints obtained using different lens models here. 

We obtain our accretio n disc const r aints using a single-epoch irn aging technique very 
similar to that presented in Bate et al. ( 2008 ) and Floyd et al. ( 20091). A similar analysis 
could (and should) be undertaken for the lightcurve technique (IKochanekl |2004| ). but is 
beyond the scope of this paper. 

As an additional benefit, this analysis offers the first opportunity for a direct com- 
parison between accretion disc parameters obtained using the lightcurve technique and 
the single-epoch imaging technique. The observational data we use, consisting of pseudo- 
broadband photometry of t he lensed quasar Q2237 +0305, was presented and analysed using 
the lightcurve technique in lEigenbrod et al.l ( 120081 ). 



- 15 - 



3.1. Method 



The single-ep och microlensing simulation technique we use here was first presented in 
Bate et al.l (|2008[ ). The goal is to constrain both the radius of the accretion disc (Tq (at a 
particular wavelength Aq), and the power- law index ( relating accretion disc radius a to 
observed wavelength A. In other words, we model the quasar accretion disc of Q2237+0305 
using the following relationship: 



a 



(9) 



where Aq is taken to be the central wavelength of the bluest filter in which observations were 
taken. 

The fu ll simulation tech nique will not be repeated here; it is described in full detail in 
Section 3 of lBate et al.l (120081 ) . There are three key differences between the analysis presented 
in that paper, and the analysis presented here. First, our Q2237+0305 magnification maps 
have a fixed microlens mass of IMq. Secondly, we use logarithmic Bayesian priors for source 
radius a, rather than the constant priors used in our previous analyses. Logarithmic priors 
ensure that the ratio of prior probability for two values of source size does not depend on 
the units chosen. This choice also allows us to compare more directly with previous work 
by other authors. And thirdly, we set the smooth matter percentage in the lensing galaxy 
at the image positions to zero. This is reasonable because the lensed images lie in the bulge 
of the lensing galaxy, where stars are expected to dominate the mass distribution. 

The result of this analysis technique is an a posteriori probability distribution for radius 
of the quasar accretion disc in the bluest filter (Tq, and power-law index ( relating accretion 
disc radius to observed wavelength. For our purposes here, we will compare the extracted 
68 per cent confidence limits on those two parameters. 



We u s ed len sing models for Q2237 -I-0305 from two sources: IWitt et al.l (Il995[ ) and 
Kochanekl (120041 ). In IWitt et al.l (119951 ) . the authors modelled four lensing systems with 
a simple family of lensing galaxy models, parameterised by the power-law index /3 of the 
lensing potential with respect to radius r. Microlensing parameters were quoted for three 
values of f3: 0.5, 1.0 (the commonly used singular isothermal sphere model) and 1.5. These 
three models allow us to test the effect of varying microlensing parameter s on ou r accretion 
disc constraints. We also use the microlensing parameters from iKochanekl (120041 ). where the 
lensing galaxy was modelled as a singular isothermal ellipsoid (SIE) with external shear. This 
allows us to directly compare results obtained usin g our single-epoch imag ing technique with 
results obtained using the lightcurve technique in lEigenbrod et al.l ( 120081 ). The parameters 
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are found in Table [31 



3.2. Observational data 



The observational data were obtained from lEigenbrod et al.l (l2008l ). who conducted 
spectroscopic monitoring of Q2237+0305 using the VLT across a three year time period. The 
data, originally FORSl spectra, were deconvolved into a broad emission line component, a 
continuum emission component, and an iron pseud o-continuum. The cont inuum emission 
was fit with a power law of the form f^, oc z/"". lEigenbrod et al.l ( 120081 ) then split the 



continuum into six wavelength bands, each with a width of 250A in the quasar rest frame, 
and integrated the continuum power-law in each band. The result is pseudo-broadband 
photometry in six w avebands, with cqntarn ination from broad emission lines and the iron 
continuum removed. lEigenbrod et al.l (120081) contai ns analysis of 39 of their 43 epochs using 
the lightcurve fitting technique of iKochanekl (120041 ). 



Here, we use two epochs of the lEigenbrod et al.l (120081 ) data, separated by a year, and 
presented in Table HI We analyse one image combination: B/A. There is no reason in 
principle why the other images could not also be analysed. Here we choose B/A both to cut 
down on simulation time, and to allow direct comparison with the lEigenbrod et al.l (l2008l ) 
results (which also only use the B/A image combination). Probability distributions for ctq 
and ( are obtained for each epoch separately, and then combined together to obtain our final 
results. 



3.3. Results and Discussion 



Results are quoted in units of the Einstein radius tiq (see Equation E]). For Q2237+0305, 
this is 9.93 x lO^^/i^n /^((M)/0.3M^)^/^cm We present 68 per c ent confidence limits on ( 
and (To for all three models in Table |5l The lEigenbrod et al.l (120081 ) results for their two cases 
(with and without a velocity prior) are provided for comparison. 

The first thing to note is that within their errors, the constraints on both (Tq and ( 
agree for all models. This suggests that at the current level of accuracy, microlensing-based 
measurements of accretion disc parameters are likely to be consistent even when different 
lensing models are used. 

However, th ere are sorae treri ds in our results which are worth discussing. If we look 
first at the three IWitt et al.l (119951 ) models, we see that as (3 is increased from 0.5 to 1.5, so 
the peak value of ( increases. It is difficult to disentangle what is causing this, however the 
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most striking variation in the magnification maps as /3 increases is an increase in the number 
of microlenses. The /3 = 1.5 magnification maps are much more densely packed with caustic 
structure than their /3 = 0.5 counterparts. 

If we comp are now the results obtained using the single-epoch imaging tec hnique (the 



2008 



Kochanekll2004j lens model) to results obtained using the lightcurve technique (the lEigenbrod et al 



results), there are a few differences. Firstly, the lightcurve technique produces smaller 
errors. This isn't surprising; we analysed two epochs of data for the single-epoch imaging 
technique, whereas Eigenbrod and collaborators analysed 39 epochs of data. Secondly, for 
the same microlensing parameters the single-epoch technique results in slightly lower (al- 
though still consistent) values for ( than the lightcurve technique. Whether this last effect 
is generically true would require more systems for comparison. 



4. Conclusion 

As a crucial step in preparing for this influx of new lensed systems that the synoptic 
survey era promises to deliver, a large-scale exploration of microlensing convergence-shear 
parameter space is warranted. In this paper, we have described how a such a parameter 
space survey could be achieved usin g a modern GPU-based supercomputing facility, and our 
experiences with the GPU-D code ( [Thompson et al.ll2010l : iBate et al.ll2010l ). 



Balancing simulation time and theoretical considerations, we suggest that the following 
is a suitable parameter set for such a survey: 

• 0.01 < K < 1.70 with increment An = 0.01 

• 0.00 < 7 < 1.70 with increment A7 = 0.01 

• 0.0 < s < 1.0 with increment As = 0.1 

• = 1000 

• A^pix = 100002 

• Map side length = 25?]o 



This results in a 170 x 171 sampling of K-7 parameter space for each s-value. The number of 
rays shot per pixel Arrays was found to be s ufficient to produce accurate magnification maps 
using the GPU-D code in lBate et al.l ( 120101 ). The pixel resolution and map size represent an 
acceptable compromise between the need to resolve the last stable orbit of a Schwarzschild 
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supermassive black hole R/scoi need for large enough maps to conduct microlensing 

simulations of quasar broad emission line regions and long light curves. 

We collated microlensing models from the hterature, and used them to suggest ways in 
which K-^ parameter space could be masked to reduce simulation time. In addition, a large 
portion of simulation time is spent in relatively few maps, where (1 — k — 7) or (1 — K + 
are close to zero. Masking out any maps with projected simulation times greater than 1 day 
using the entire facility enables swifter completion of the survey, at the risk of missing some 
interesting areas of parameter space. With these restrictions in place, we estimate a total 
simulation time of ~ 414 days utilising the entirety of a 100 teraflop/s facility (gSTARioo), 
or ~ 167 days on a 600 terafiop/s facilty (gSTAReoo)- Realistic scheduling requirements for 
a high-demand supercomputer will increase the total wall-time required, so a staged strategy 
starting with parameters for the currently known microlensed quasars, would support early 
science usage. 

Finally, we have provided an example of one possible application for a large microlensing 
parameter survey. We used four different microlensing parameter sets to take a preliminary 
look at the systematic effects of varying lens models on accretion disc constraints in the lensed 
quasar Q2237-I-0305. Although all of the lens models produce results that are consistent 
within their errors, we note that there are suggestions that models with denser microlensing 
environments lead to shallower accretion disc temperature profile predictions. 

This research was supported under Australian Research Council's Discovery Projects 
funding scheme (project number DP0665574). We are grateful to Alex Thompson for his 
work on the GPU-D code. Our referee provided very helpful suggestions, that improved the 
quality of the work. 

Facilities: VLT:Kueyen (FORSl). 
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5 10 
r?o (10^^ No"^^^ [<M>/0.3Mq]^/^ cm) 

Fig. 1. — The Einstein radius projected on to the source plane rjQ for each of the 59 
lensing systems with both source and lens redshifts on the CASTLES Survey webpage 
|http://cfa-www.harvard.edu/glensdata/). The Einstein radius is calculated for a O.SMq 
lens. The solid line indicates the mean r]o, and the dashed lines the standard deviation. 
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7 



Fig. 2. — A contour plot showing simulation time for a magnification map with a given 
combination of n and 7. The colour gradient ranges from hours (black) to 24 hours 
(white). Overlaid are the existing models taken from the literature (Tabled]). Times were 
calculated using Equation IH assuming maps with a side-length of 2brjQ, N^i^ = 10, 000^ and 
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Fig. 3. — A contour plot showing simulation time for a magnification map with a given 
combination of k and 7, with maps that require simulation time greater than 1 gSTARioo day 
masked out. The colour gradient ranges from hours (black) to 24 hours (white). Overlaid 
are the existing models taken from the literature (Tabled]). Times were calculated using 
Equation [HI assuming maps with a side-length of 2brjQ, N^i^ = 10, 000^ and A^, 



' rays 



1000. 
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Table 1. Known microlens models and references. 



System 


-^models 


Ref. 


HE 0047-1756 


1 


(4) 


QJ 0158-4325 


10 


(7) 


MG 0414+0534 


3 


(13) 


HE 0435-1223 


1 


(4) 


HE 0512-3329 


1 


(4) 


SDSS 0806+2006 


1 


(4) 


SDSS 0909+532 


2 


(4), (5) 


SDSS J0924+0219 


13 


(2), (4), (6) 


FBQ 0951+2635 


1 


(4) 


QSO 0957+561 


1 


(4) 


SDSS JlOOl+5027 


1 


(4) 


SDSS J1004+4112 


1 


(4) 


QSO 1017-207 


1 


(4) 


HE 1 1 04-1 805 


1 1 

X X 


(A) (7) 


PG 1115+080 


5 


(4), (8), (13) 


RX J1131-1231 


11 


(1), (4) 


SDSS J1206+4332 


1 


(4) 


SDSS J1353+1138 


1 


(4) 


H1413+117 


4 


(4), (13) 


B J1422+231 


1 


(4) 


SBS 1520+530 


1 


(4) 


WFI J2033-4723 


1 


(4) 


Q2237+0305 


12 


(3), (9), (10), (11), (12), (13) 



References. — ('. 


) 


Dai et al. 2010: (2) 


Keeton et al. 


2006 


(3) 


K^ochanek 




2004: (A) Mediavilla et al. 


2009: 


(5) 


Mediavil 


a et al. 




2011: f6) iMorean et al. 


booe 


i 


Morgan et al. 


2008 


(8) 


Poolev et al..,2009£ (9) Rix et al 


1992: 


(10) 


Schmidt et al 




1998: (11) Schneider et al. 


1988 


(12) 
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Wambsganss fc Paczvnskilll994l : flSl lWitt et alJll995 



Table 2. Simulation strategies. 



Strategy 


Side length (rjo) 


Side length (pixels) 


Time (gSTARioo) 


Time (gSTAReoo) 


Notes 


Complete 


50 


16500 


115 years 


19 years 


full parameter space simulation 


Realistic 


25 


10000 


414 days 


167 days 




Stage 1 


25 


10000 


122 days 


40 days 


one map at each Table [T] model position 


Stage 2 


25 


10000 


177 days 


59 days 


maps at Table [T] model positions ±0.01 


Stage 3 


25 


10000 


332 days 


140 days 


maps at Table [T] model positions ±0.05 


Stage 4 


25 


10000 


392 days 


157 days 


maps at Table [T] model positions ±0.1 



Note. — Parameter space is tiled with a resolution of 0.01 in both re and 7, and and 0.1 in s. All times are estimated, and assume 100% 
utilisation of gSTAR. The 'Complete' strategy assumes no time constraints whatsoever, and is obviously impractical. All subsequent strategies 
exclude any parameter combinations that require greater than 1 device day per map to simulate, and mask out regions of parameter space 
containing no models (see Figure [3]l . Stages 1-4 represent progressive filling of parameter space towards the Realistic strategy. The listed 
simulation times are cumulative - for example, Stage 2 takes 177 gSTARioo days to complete, but at the end of that time period Stage 1 will 
also have been simulated. 
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Table 3. Microlensing parameters for Q2237+0305. 



T\ T 11 

Model 


Image 


K 


7 


Magnification fi 


Ref. 


(3 = 0.5 


A 


0.228 


0.603 


4.3 


(1) 




B 


0.223 


0.582 


3.8 






C 


0.289 


0.964 


-2.4 






D 


0.261 


0.880 


-4.4 




/3 = 1.0 


A 


0.471 


0.416 


9.4 


(1) 




B 


0.464 


0.406 


8.1 






C 


0.552 


0.617 


-5.6 






D 


0.515 


0.579 


-10.0 




13 = 1.5 


A 


0.728 


0.215 


36.5 


(1) 




B 


0.723 


0.212 


31.3 






C 


0.789 


0.296 


-23.3 






D 


0.761 


0.286 


-40.7 




SIE+7 


A 


0.394 


0.395 


4.7 


(2) 




B 


0.375 


0.390 


4.2 






C 


0.743 


0.733 


-2.1 






D 


0.635 


0.623 


-3.9 





References. — fll lWitt et al.lll995l : {2) iKochaneM 12004 
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Table 4. Two epochs of observational B/A and C /A flux ratios for Q2237+0305. 



Band Emitted Ac (A) Observed B/A C/A Date 



1 1625 ± 125 4379 ±337 0.52 ± 0.02 0.31 ± 0.03 2005 November 11 

2 1875 ± 125 5053 ±337 0.51 ±0.02 0.33 ± 0.03 2005 November 11 

3 2125 ± 125 5727 ±337 0.50 ±0.01 0.34 ± 0.03 2005 November 11 

4 2375 ± 125 6401 ± 337 0.49 ± 0.01 0.35 ± 0.03 2005 November 11 

5 2625 ± 125 7074 ± 337 0.48 ±0.01 0.37 ±0.03 2005 November 11 

6 2875 ± 125 7748 ± 337 0.47 ±0.01 0.37 ±0.02 2005 November 11 

1 1625 ± 125 4379 ± 337 0.33 ± 0.02 0.28 ± 0.02 2006 November 10 

2 1875 ± 125 5053 ± 337 0.34 ± 0.02 0.30 ± 0.02 2006 November 10 

3 2125 ± 125 5725 ± 337 0.35 ± 0.01 0.31 ± 0.02 2006 November 10 

4 2375 ± 125 6401 ± 337 0.36 ± 0.02 0.33 ± 0.02 2006 November 10 

5 2625 ± 125 7074 ± 337 0.37 ± 0.02 0.34 ± 0.02 2006 November 10 

6 2875 ± 125 7748 ± 337 0.37 ± 0.01 0.35 ± 0.02 2006 November 10 



Note. — Data obtained from lEigenbrod et al.l (120081 ). Photometry was extracted from 
spectra obtained with the FORSl spectrograph on the Very Large Telescope (VLT) at the 
European Southern Observatory (ESO). Filter wavelengths are i n the rest frame of the 
source quasar, located at a redshift of Zs = 1.695. Following the lEigenbrod et al.l ( 120081 ) 
numbering, the 2005 November 11 dataset is epoch number 17, and the 2006 November 
10 dataset is epoch number 28. 
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Table 5. Accretion disc constraints for three Q2237+0305 models. 



Model Radius ctq Power-law index ( 



f3 = 0.5 




< 0.15770 


0.79^°;P 


/3 = 1.0 




< O.lbrjo 


0.88t°:f7 


13 = 1.5 




< o.ier^o 


1 40+0.83 


Kochanek r2004) 


model 


< 0.21r/o 


' ' -0.34 


Eieenbrod et al. 


^2008) velocity prior 




1.2 ±0.3 


Eieenbrod et al. 


^20081 no velocity prior 




1.1 ±0.3 



Note. — Constraints on radius of the accretion disc in Q2237±0305 at 
^obs = 4379A, and power-law index relating accretion disc radius and ob- 
served wavelength of emission. C onstraints were obtaine d using the B/A 
image combination. Results from lEigenbrod et al.l (120081) . which were ob- 
tained using the lightcurve technique of iKochaneld (l2004l ). are provided for 
comparison. 



